Quantum simulation of cosmic inflation in two-component Bose-Einstein condensates 
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Generalizing the one-component case, we demonstrate that the propagation of sound waves in two- 
component Bose-Einstein condensates can also be described in terms of effective sonic geometries 
under appropriate conditions. In comparison with the one-component case, the two-component setup 
offers more flexibility and several advantages. In view of these advantages, we propose an experiment 
in which the evolution of the inflaton field, and thereby the generation of density fluctuations in the 
very early stages of our universe during inflation, can be simulated, realizing a quantum simulation 
via analogue gravity models. 
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I. INTRODUCTION 



Within our present standard model of cosmology, basi- 
cally all inhomogeneities ~ including the seeds for the for- 
mation of structures such as our galaxy - originate from 
quantum fluctuations of a single scalar field, the infla- 
ton [EH- This (postulated) field drives inflation, which 
is a stage of very rapid expansion in the earliest evolu- 
tionary phase of our universe |3]- Tracing the inflaton 
fluctuations back in time and thereby undoing the red- 
shift induced by the cosmic expansion, the anisotropies 
of the cosmic microwave background we observe today 
correspond to extremely short wavelengths during infla- 
tion. As a result, the fluctuations of the cosmic mi- 
crowave background probe ultra-high (e.g., Planckian) 
energy scales - which are experimentally inaccessible 
with the present-day and near-future available technol- 
ogy of, for example, particle accelerators. At such ultra- 
high energies, quantum effects of gravity are expected to 
become important - but the underlying physical theory 
for the description of these effects is not known yet. Con- 
sequently, high-precision measurements of the cosmic mi- 
crowave background might give us some insight into new 
physics beyond well-established theories (in particular, 
beyond the standard model of particle physics). 

In order to detect signatures of the new physics in the 
anisotropies of the cosmic microwave background, one 
has to investigate which kind of higher-order corrections 
and correlations could potentially be induced by devia- 
tions from the known laws of physics occurring at ultra- 
high energies. One way to achieve this aim is to consider 
analogous systems, based on laboratory physics, which 
reproduce major features of the inflaton field, and which 
can therefore be used to simulate the generation of (quan- 
tum) fiuctuations and further interesting effects - theo- 
retically as well as experimentally. This line of approach 
has come to be known under the term analogue grav- 
ity/cosmology, see, e.g., |^. The consideration of these 
analogues leads to a better understanding of the system 
to be simulated, in particular, regarding the possible im- 
pact of high-energy degrees of freedom. As cosmology is 



essentially a purely observational branch of science, be- 
cause obviously we cannot do experiments on the real 
system, i.e., the universe, the analogues also allow an 
experimental verification of so far only theoretically pre- 
dicted effects. Basically the same motivations underlie 
the idea of quantum black hole analogues ("dumb holes" 
[a, la I3)> since the Hawking radiation of "real" grav- 
ity black holes has its origin in (trans-)Planckian modes: 
The quanta emitted necessarily come from regions very 
close to the black hole horizon, and experience a large 
redshift when finally detected far away from the black 
hole 5, 8]. 

A scalar field (such as the infiaton) within the 
curved space-time structure of an expanding universe 
can be simulated by propagating sound waves in single- 
component Bose-Einstein condensates (BEC) Q; pre- 
vious considerations on such effective acoustic geome- 
tries in single-component BECs can be found in Refs. 
HllSlillillil- The advantage of BECs lies in the fact 
that the corresponding parameters (such as local density 
and speed of sound) can be controlled with atomic pre- 
cision experimentally, and that the underlying physics is 
well-understood on all energy scales in particular in the 
dilute case, for which the Gross-Pitaevskii equation pro- 
vides a rather accurate description of the order parameter 
dynamics. There are basically two possibilities for simu- 
lating an expanding universe within a BEC: changing the 
interparticle coupling or expanding the condensate (or a 
combination of both) 11, 12, 13] . However, both meth- 
ods come with problems. Firstly, a controlled expan- 
sion of the BEC cloud requires a specific time-dependent 
trap; and furthermore, the density of the cloud rapidly 
decreases during the expansion, leaving only a short time 
to do the experiment. Secondly, in order to change the 
inter-particle coupling drastically via a Feshbach reso- 
nance, by rapidly sweeping a time-dependent external 
magnetic field, one has to go very close to the reso- 
nance; one then encounters the problem that the cou- 
pling constant effectively acquires an imaginary part due 
to molecule formation in three-body recombination pro- 
cesses [lj,[i3, which spoils the desired effect. Thirdly, a 
time-dependent interparticle coupling would also induce 



a variation of the density and the size of the BEC cloud 
- unless the trap is changed accordingly, which again is 
difficult (see first point). Finally, one has to be able to 
measure the generated fluctuations in order to simulate 
the inflaton field experimentally. 

In the present article, we investigate whether it is 
possible to overcome some of these problems in two- 
component Bose-Einstein condensates, which are readily 
experimentally available; the various components can be 
realized by trapping different hyperfine ground states of 
the same atom U^ [13, [lE H^ ■ We start by describing in 
scctionlnlhow the collective equations of motion for small 
fluctuations (i.e., sound waves) in the two-component gas 
are obtained. In section ITTTl we show that it is possible 
to map these equations of motion onto two effective sonic 
metrics (bi-metricity) under certain conditions. The ana- 
logue of the Planck scale, where the concept of the ef- 
fective geometry breaks down, is discussed in section HVI 
Section^is devoted to the simulation of the de Sitter ge- 
ometry - i.e., to constructing an analogue for the inflaton 
field during inflation. We discuss possible experimental 
realizations and measurement prescriptions to realize the 
desired inflation quantum simulation in section IVll The 
advantages and drawbacks of the described method and 
further aspects are summarized in section IVIII 

We note that the inflaton has been mentioned in 
the context of Bose-Einstein condensates previously [23| ■ 
However, the discussion there has been rather qualita- 
tive. In particular, the inflaton mode has not been re- 
lated there to any effective space-time metric of cosmo- 
logical character (for example the de Sitter metric), in 
which its propagation ought to take place. The latter is 
necessary to appropriately describe the freezing process 
of the quantum fluctuations and the related concept of a 
horizon, both of which we shall investigate for de Sitter 
space-time in what follows. 



II. EQUATIONS OF MOTION 

In terms of the Madelung representation for the order 
parameter components 7/'a(T,t) 



^pa{r,t) = y/gjy^ eyip{iSair,t)} , 



(1) 



with the density Qair, t) and the phase (eikonal) Sa{r, t), 
the Lagrangian density of a dilute two-component Bose- 
Einstein condensate reads (we put h ~ 1) 
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Here, the two masses of the atoms are nia {a = 1,2), 
the one-particle trapping potentials (which are generally 
different) are given by V^(r, i), and the (symmetric) two- 
particle interaction coupling matrix is denoted gab- 



Linearizing around a given, stationary background so- 
lution e°(r) and S^{r) with v^ = VS^/ma and neglect- 
ing the quantum pressure terms ex (V^ftT)^ - which 
amounts to the local density (Thomas-Fermi) approxi- 
mation in one-component Bose-Einstein condensates - 
leads to the second-order effective action 
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Varying the above action with respect to Sga yields two 
Bernoulli-type equations for the fluctuations 



DaSSa + 22 9abSgb = , 



(4) 



where the co-moving derivative is deflned to be 

DaSSa = dtSSa + V^ ■ V5Sa . (5) 

(2) 

Using the above equation to eliminate 5ga from C^g , we 
obtain a phases-only effective Lagrangian of the form 



ab 



2ma 



The general wave equations for the phase fluctuations 
dSa then take the form 
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J2DA9-b'DbSS,)-^V'SSa = 



(7) 



These wave equations, in the general case, do not yet 
have the pseudo-Lorentz invariance required to obtain 
effective space-time metrics of Lorentzian signature. 



III. EFFECTIVE GEOMETRY 

So far, we have involved no specific assumptions about 
the background densities and velocities as well as the 
constituent masses. We now come to discuss a simple 
case in which an effective metric descri ptio n in terms of 
the Painleve-GuUstrand-Lemaitre type J2l| is viable. 

As demonstrated in Ref. |22|, the introduction of ef- 
fective geometries for multiple interacting fields is more 
involved than the single-field case - where rather general 
assumptions ensure the existence of an effective metric 
for the propagation of perturbations. In the case of mul- 
tiple interacting fields, there are the following three main 
possibilities: a) Owing to a lack of symmetry one cannot 
introduce a metric at all ("pre-geometry"); b) the per- 
turbations effectively decouple and can be described by 
multiple metrics; or c) all the metrics coincide and there 
is one unique metric (reproducing the principle of equiv- 
alence). There is also the (fourth) possibility that the 



propagation of perturbations is not equivalent to scalar 
(i.e., spin-zero) fields in curved space-times, but to fields 
with higher, non-zero spin instead (e.g., Dirac |2.1| or vec- 
tor 24] fields). 

Phonons in arbitrary two-component condensates cor- 
respond to case a) in general - even though one might 
diagonalizc the dispersion relation, the full equations of 
motion do not allow the introduction of an effective ge- 
ometry (in the most general situation) . In order to arrive 
at an effective metric, certain requirements on the back- 
ground solution are necessary. Let us assume that the 
parameters of the background solution satisfy the follow- 
ing conditions 



and the eigenvectors read 



SSi ± SS2 

V2 



(12) 



If, in addition, .gdiag ~ <?offj which can be fulfilled to 
a high degree of accuracy (on the level of three percent) 
between different hyperfine species in *TR,b (Ref. [l3|) as 
well as in ^■^a (Ref. U9i|), we have one "hard" density 
mode (f)+ , and one "soft" spin mode 0_ . This separation 
of energy scales occurs close to the point of spatial phase 
separation of the two components, due to the increased 
interspecies repulsion g^ff 
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eigenvectors (gafc) = 


= const. 





(8) 



The diagonalization of the (real and symmetric) coupling 
matrix gab leads to eigenvalues g± , given by 



9± 



911 + 922 



± 



.911 - .922 



.9?2 



(9) 



Note that ~ in contrast to the eigenvectors of the ma- 
trix gab - its eigenvalues g± are not required to be con- 
stant. By virtue of the assumptions (jHJ, the Lagrangian 
in Eq. © can be diagonalized 



C, 



■off 



1 



= ^\/-&±S±df,<j)±d,,(l)±, 



(10) 
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where (j)± denote the projections of the phase fluctuations 
dSa onto the (constant) eigenvectors of gab and a sum- 
mation convention over indices /i, v and ± is implied. In 
this (highly symmetric) case, we therefore obtain two in- 
dependently propagating, i.e., decoupled modes (f>±, which 
feel effective space-time metrics of the conventional (co- 
variant) Painleve-GuUstrand-Lemaitre form [23 



:ii) 



where the two sound velocities are c± — y^g±go/m and 
1 represents the unit matrix. As long as these two sound 
velocities do not coincide c+ 7^ c_ , the system under con- 
sideration corresponds to the bi-metric case b) discussed 
at the beginning of this Section. 

The assumption in JSJ that the eigenvectors of gab be 
constant can be satisfied, if gafc itself is constant, or if it is 
sufficiently symmetric. We shall assume gu — (722 in the 
following, because it allows both for a bi-metric approach 
and an implementation of time-dependent g± = g±{t). 
In this situation of gu == g22 = <7diag, the eigenvalues are 
simply given by g± = gdiag ±5off, where 512 = 321 = ffoff, 



IV. DISPERSION RELATION AND THE 
PLANCK SCALE 

So far we discussed the effective geometry for low- 
energy excitations (sound waves). The full dispersion 
relation, without the Thomas- Fermi approximation (i.e., 
not neglecting the quantum pressure terms) can be ob- 
tained via the JWKB approximation, which amounts to 
the geometrical optics limit of quasiparticle propagation. 

Combining the (linearized) equation of continuity 



t{Lo + v'i-k) 6ga = — k^SSa , 



(13) 



and the Bernoulli-type equation Q augmented with the 
quantum pressure term on the right-hand side, 



t{uj + v'i-k) SSa + ^ gabSgb = -j-—n 

'illiaiJa 



Sga, (14) 



gives us two Bogoliubov dispersion relations of the usual 
type. Using the assumptions ((SJl, we obtain 



m 2 fc 

— (cj± +v-k) + 

go 4:171 go 



= 9±k' 



(15) 



The deviation from the linear (rest frame) dispersion 
w^ oc k^ occurs at the two healing lengths 



e. 
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Am gag ± 



(16) 



Accordingly, the analogues of the Planck length, the two 
coherence lengths S,± , scale in the same way as the inverse 
spin and density mode velocities. 



e^ci 
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(17) 



V. 



DE SITTER SPACE-TIME AND ANALOGUE 

INFLATON 



The extreme dependence of the sound velocity of the 
spin mode on the coupling matrix gab can be profitably 



used to simulate a rapidly expanding universe via small 
temporal changes gab — dabit)- A "spin horizon" for 
the spin mode should be easier to realize experimentally 
than the sound horizon of a one-component BEC, in view 
of the possibility to manipulate the "spin" velocity such 
that it closely approaches zero. 

The line elements for a background at rest, v = 0, read 



dsl 



go c± dt 



1 



dr' 



(18) 



From now on we focus on one particular mode, the spin 
mode, drop the subscripts ± in most of the following 
formulas, and furthermore set to = 1 for convenience. We 
suppose the background density go to remain essentially 
constant during rapid variations of g- (which is possible 
for Va = Vb, cf. section |VII|) . 

Assuming a time dependence of the propagation veloc- 
ities/coupling constants 



Co 
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(19) 



with H being the condensed-matter analogue of the Hub- 
ble parameter in cosmology, we obtain the de Sitter met- 
ric 



ds^ — goCo dr^ — 



^Ht 



Co 



dr^ 



(20) 



with a transformed de Sitter time-coordinate r = 
H~^\a{Ht), representing proper time (not equal to the 
laboratory time), and the prefactor goCo = const. 
The corresponding Klein-Fock-Gordon equation 



Ud)^ 
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■d,, 



ra.</))=o. 



(21) 



with g_ denoting the (negative) determinant of the met- 
ric g~ and g^ its inverse, assumes in this case the simple 
form (in three spatial dimensions) 



dr^ dr 



-2Ht 



[co^Y U - 



(22) 



After a spatial mode expansion into plane waves (tak- 
ing into account isotropy and homogeneity), each mode 
behaves as a damped harmonic oscillator with a time- 
dependent potential 



• horizon-crossing: At some point in time, the 
monotonously decreasing term e~^^'^ , correspond- 
ing to the expansion of the universe, becomes small 
enough and the damping term starts to play a role. 
Since the length scale cq/H corresponds to the size 
of the particle horizon in an expanding (de Sitter) 
universe, this is the point where the modes cross the 
horizon and hence do not oscillate freely anymore. 

• freezing: For late times, we have e~^^'^[cofe]^ ^ 
H^, and hence the potential term can be neglected. 
This corresponds to a strongly over-damped oscil- 
lator and thus the modes do effectively not evolve 
anymore. 

The same three essential stages undergoes the inflaton 
field during the epoch of inflation in the very early uni- 
verse, according to our (present) standard model of cos- 
mology j2|. Since the Hubble parameter during inflation 
is expected to be significantly smaller than the Planck 
scale - otherwise the (semi-classical) notion of a (quan- 
tum) field within a curved space-time would not apply - 
the modes leave the Planck regime (are "born" ) nearly 
freely oscillating. Later on, with increasing r, they cross 
the horizon and freeze. At the end of inflation, the frozen 
quantum fluctuations of the inflaton field are transferred 
into (classical) density fluctuations (decoherence of the 
quantum state due to interaction with other degrees of 
freedom) , which are in turn supposed to be the seeds for 
structure formation in our universe represented, e.g., by 
our galaxy. 



VI. MEASUREMENT OF FLUCTUATIONS 

In order to discuss the quantum state of the fluctua- 
tions (f), it is convenient to introduce yet another time- 
coordinate, i.e., the conformal time rj — —e^^'^/H = 
— l/(7J^i), in terms of which the de Sitter metric H20() 
can be cast into the conformally flat form 



ds^ ^ 



goCQ 



drf 2 dr" 



(24) 



Expanding the phase operator (j) into plane-wave solu- 
tions of the Klein-Fock-Gordon equation H21() in confor- 
mal time 
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(25) 
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-2Ht 



cofc] Ufc = 0. (23) 



The evolution induced by the above equation of motion 
can roughly be split up into three regimes: 

• oscillation: For very early times, we have 
g-2ffr|-^^^j2 ^ ^2. j^gj^Qg ^jj^g damping term can 



we have, for an arbitrarily chosen quantization volume 
V, the analytical expression 
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(26) 



We introduced creation and annihilation operators a^ 
be neglected and the modes oscillate almost freely. and a^, where the dk annihilate the "adiabatic" vacuum 



state |0)^^ (see, e.g., |2g) 



afe|0) 



ad 



0. 



(27) 



For early times (in the oscillating regime rj | — oo), the 
adiabatic quantum vacuum state coincides to zeroth or- 
der with the instantaneous ground state and is therefore 
a natural candidate for the vacuum state. As indicated 
above, after horizon crossing, the fluctuations are frozen 
at late times ifoo-^^ryfO. In this regime, we obtain 
from (|26l) the expectation value 



2Fc3fc3 



m'- 



(28) 



We display the evolution of a (j)^ component of H2()l) and 
the approach to the above "frozen" value in Fig. Q. 

In complete analogy to the cosmic inflaton field, the 
expression (I28II generates a scale-invariant spectrum: In 
view of the d^fc-integration, the corresponding spatial 
correlation function is (within the region of validity of 
the above approximations) invariant under re-scaling 
(<^_(r)<^_(r')) c, (0_(Ar)<^_(Ar')), see, e.g., ^. 



lm[0fc(t)] 




FIG. 1: Freezing process for the modes (f>k according to ex- 
pression 12611 , with the absolute value of the frozen <^^ given 
by I28I I. displayed on a logarithmic lab time scale t. 



A. Phase-phase correlations 

Let us investigate the impact of the frozen quantum 
fluctuations on the (spatially) Fourier-transformed two- 
point phase-phase correlation function defined as 

C^{k,t)^ f (fr cxp{ik-r}{$^{0,t)$^{r,t)). (29) 

Inserting Eq. (|28|l we obtain (after freezing) 

H'go 



C^{k) 



2c3fc3 



(30) 



It appears that one could generate an arbitrarily strong 
effect by increasing the Hubble parameter H , but this 



is a fallacy. The validity of the approximations made 
requires that the conditions c^k 3> -ff S> cgnai^ on the 
frequency scales hold (with Cfinai denoting the final value 
of the time-dependent speed of sound), as well as that 
the long wavelength limit k£^-{t) ^ 1 is fulfilled. In or- 
der to estimate the maximally obtainable effect, we insert 
iJ ^ cofc as well as fc ~ k^lt^'^^ '^itl^ f^R^tiS'^^ denoting the 
inverse of the final healing length, i.e., the final Planck 
scale. The maximum effect of the dimensionless quan- 
tity {Ak)^C^{k), i.e., the correlation within some wave- 
number interval (Ak)^, then reads 



(Afc)3c^^'^(fc) 
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T^Planck 
'Cfinai 



(31) 



Since the concept of the effective geometry only applies 
for k <C fcfina'i"^'*^' ^^^ maximum impact of the frozen 
quantum fluctuations on the two-point phase-phase cor- 
relation function is limited by the product of the ra- 
tio vyoTfffinai 3> 1 and the final diluteness parameter 
V^PoSflnai ~ which is bound to be a quite small quan- 
tity; for the spin mode of relevance here it is of order 
10-2... 10-4. 



B. Density-density correlations 

In contrast to phase-phase correlations, which can 
only be measured by time-of-flight measurements |29j, 
density-density correlations can be obtained directly, i.e., 
in situ, for example via (state-selective) absorption (such 
in situ measurements are generally more difficult to be 
carried out, though [23 ). 

The spatially Fourier-transformed two-point density- 
density correlation function is defined, in a manner anal- 
ogous to its phase-phase counterpart in H29|l . to be 



C,{k,t) 



/ (fr exp{ik ■ r}{dg^{0,t)5g^ 



(r,i)). (32) 



Note that the expectation values in Eqs. H29I) and 1)32(1 
are realized via statistical averages in any measurement 
since the quantum state is not an eigenstate. 

To calculate the density fiuctuations, we have to derive 
the time-dependence of the phase fluctuations beyond the 
zeroth-order (frozen) part (f)^. The approximate solution 
of Eq. 1(23(1 respectively Eq. 1(25(1 for frozen inflaton modes, 
transforming back to laboratory time, reads 



Mt)^<t>li^ + 



2HH^ 



-'i^ 



(33) 



Since g ex l/f*, the density fluctuations increase linearly 
in time, cf. Eq. Q): 
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90 



(34) 



The overdot denotes a partial derivative with respect to 
lab time. The sub-leading term O (i°) ensures the valid- 
ity of the canonical commutation relation for hydrody- 
namical density and phase operators, [Sg{t, r), (f>{t, r')] = 
iS{r — r'). 

Using the above relations (jM|l and (^5)) . and inserting 
into (|23) yields 



C,{k,t) 



250 



(35) 



In addition to the assumptions already discussed in the 
case of phase-phase correlations (cofc S> i? ^ cgnaifc and 
A:^_(t) ^ 1), the total time duration of the sweep is 
limited by the maximally possible variation of the cou- 
pling t^jjg^j = co/{cfini,iH^), according to Eq. ^. Again 
we estimate the maximally obtainable effect by setting 
H ^ cok, k ~ fcPltr^ and i^„^i = co/(cfi„aii?'). In- 
serting into Eq. (|35l) . the maximally obtainable relative 
change induced in the two-point density-density correla- 
tion function in the wave-number interval Ak reads 



(Afc)3C"-'^(fc,ifl„al) 
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,(36) 
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where the last relation also uses, from 

V^fi'Offfinal, that H fa Ctinal\/Q090 — CfinalCo- 

The density-density correlations are increasing with 
time. Howe ver, the maximal relative change is still by 
a factor of \/<?o/ff final 3> 1 smaller than in the case of the 
phase-phase correlations, Eq. 1)31(1 . For density-density as 
well as phase-phase quantum correlations, the final out- 
come is rather small owing to the intrinsic smallness of 
the final diluteness parameter \/qo9^^i- 



C. Amplification in unstable regime 

It became evident in the preceding considerations that 
the impact of the frozen vacuum fluctuations is in princi- 
ple observable, but rather weak. Therefore, it is desirable 
to find some mechanism for the amplification of the fluc- 
tuations before measuring them. As one possibility, we 
propose changing the sign of g- after the inflation phase, 
i.e., switching to the unstable regime, for a short period 
of time. If this change in g- occurs much faster than the 
(frozen) dynamics of the field (p- , we may apply the sud- 
den approximation, assuming a step-function like change 
of g- from a small and positive g"^ to a small negative 
value 5°"'. 

Within the sudden approximation, the relation of the 
behavior of the field (/)_ just before (0'i\ (/>'") and after 



can be obtained from the equations of motion 



r(r,t) = <^°"'(r,t), 
1 



i^'-(^><) = ^^-'ir,t) 



9- 



5- 



(37) 



The latter equality is equivalent to 6g^^{r,t) — 
(5p°"*(r,i). Assuming a constant g™* < after the tran- 
sition to the unstable region, the subsequent time evo- 
lution is simply given by the superposition of the two 
independent solutions of imaginary spin mode frequency 



i>r 



(t) = Ak exp 



-Bk exp 




(38) 



Fluctuations with large wavenumbers (but still small 
compared to the Planck scale) grow faster. Note that this 
behavior is opposite to that of the gravitational (Jeans) 
instability in the early universe, where small wave- 
numbers experience a stronger amplification ~ larger 
structures collapse faster due to the gravitational at- 
traction than inhomogeneities on smaller length scales 
27, 281. 



If we had (/-^"V^^' 



eo|g°"'|fe^/m, the factor 



iout J^OUt \ 



the rapid change of g- from g^ > to g™* < 



Ak and hence the growing component of the solution 
would exactly vanish. However, during the phase of freez- 
ing, the time-derivative (j)k decreases and is finally sup- 
pressed by a factor of yoinai/co- Hence the exponen- 
tially increasing part Ak does contribute in a roughly 
equal amount, Ak f^ Bk- As a result, we obtain a drastic 
amplification of the phase (and also density) fiuctuations 
until the non-linear regime (that is, phase separation 17j) 
is reached, where the fiuctuations are not small compared 
to the order parameter itself anymore. Consequently, 
one would expect the fluctuations to be measurable - 
with state-selective absorption imaging, for example - 
provided that all other perturbations are small enough 
(see the next section). 



VII. CONCLUSIONS 

In summary, it is possible to simulate the behavior of 
the infiaton field - more accurately, its quantum fluctua- 
tions - within two-component Bose-Einstein condensates. 
If we decrease (?_ with time according to Eq. (|19|) . the 
phase fluctuations (/)_ (i.e., the spin mode) experience an 
effective de Sitter metric, and hence undergo the three 
stages of free oscillation, horizon crossing, and freezing 
explained in section The frozen initial quantum fluc- 
tuations can be measured after an amplification phase, 
in which the coupling matrix is tuned to 5_ < 0. 

The difference of the spectra (in fc-space) of the frozen 
initial quantum fluctuations (with an approximately scale 



invariant spectrum ex 1/fc'^) and the usual quantum fluc- 
tuations of a condensate at rest (ex 1/A;), provides one 
possibility to distinguish the signal from other effects. As 
another option for a consistency check, one could keep the 
coupling g- > constant for a given time duration after 
the de Sitter phase (which has a decreasing 5^ ex t^'^ > 0) 
and before switching to the unstable regime g- < 0. This 
way, one would allow the modes to freely oscillate again 
before amplification, and thereby suppress those modes 
which are in the second or fourth quarter of their pe- 
riod (i.e., have a decreasing amplitude) at the moment 
of switching to the unstable regime. This intermediate 
phase of small constant f;_ > would be roughly anal- 
ogous to the epoch after inflation, and the suppression 
mechanism is similar to the suppression of modes in the 
cosmic microwave background, leading to the peaks and 
valleys in its power spectrum p30| . 



A. Advantages 



ters - which can be resolved with optical methods (e.g., 
state-selective resonant absorption imaging). 

Furthermore, the measurability of the quantum fluc- 
tuations requires the absence of any other fluctuations 
(noise) which may swamp the signal. Assuming an initial 
speed of sound of about 1 mm/sec, a mode with a wave- 
length of a few fiia corresponds to an initial frequency 
somewhere below 1 kHz. Hence one should avoid temper- 
atures above 100 nK, for which the thermal fluctuations 
exceed the (to be measured) quantum fluctuations (as- 
suming no thermalization during the experiment). For 
the parameters specified above, the Hubble parameter 
determining the timescale for changing g- would also be 
below 1 kHz. This should represent no problem as it is 
possible to change the coupling constants (i.e., the exter- 
nal magnetic field) on the time-scale of some microsec- 
onds. 

Although quite challenging, the proposed experiment 
should thus in principle be feasible with present-day tech- 
nology. 



In comparison with the one-component case, the re- 
alization of an effective geometry simulating the expan- 
sion of the universe in a two-component Bose-Einstein 
condensate has several advantages: Firstly, in order to 
generate large relative variations of the spin coupling f/_ , 
it is not necessary to go near the Feshbach resonance, 
and therefore the problems associated with molecule for- 
mation (inducing an imaginary part in the coupling con- 
stant) mentioned in the introduction are not relevant. 
Secondly, it is possible to change g- while sustaining a 
constant background density: Assuming Va = VJ, and 
g+ = const, we obtain v^ — Q and g^ = 0, g^ = as 
a valid solution even for time-dependent g~{t). Thirdly, 
it is possible to amplify the fluctuations in the unstable 
regime g^ without necessarily destroying the condensate 
by phase separation. 



B. Experimental parameters 

In order to discuss the experimental feasibility of the 
proposed quantum simulation of the inflaton, it is neces- 
sary to provide an estimate for the experimental param- 
eters involved. If we assume that we may achieve essen- 
tially equal initial healing lengths for the spin and the 
density mode ^i!? « ^™ = 0(100 nm), we may simulate 
an expansion of the universe corresponding to more than 
one e-folding (i.e., HAt > 1), arriving at the final healing 
lengths C"* = 0(1 ^m) and ^°"* = ^^ = 0(100 nm). Al- 
though one e-folding is tiny compared to the expansion 
during inflation (which involves many orders of magni- 
tude), one can still reproduce generic features (such as 
freezing). Since the concept of the effective geometry 
only applies for wavelengths large compared to the heal- 
ing length (the analogue of the Planck scale), the charac- 
teristic size of the interesting structures (and hence the 
size of the condensate cloud) should be several microme- 



C. Planckian problem 

Apart from simulating the inflaton fleld experimentally 
(quantum simulation via analogues), one of the long-term 
motivations of the present work was to be able to ulti- 
mately investigate the impact of the behavior at the effec- 
tive Planck scale on the frozen fluctuations (as explained 
in the introduction). For dilute Bose-Einstein conden- 
sates, the behavior for small wavelengths (i.e., below the 
healing length) is in some sense trivial - the excitations 
are free particles of mass m (the constituent atoms), with 
Lu — k /{2m). More importantly, though, the healing 
length increases with time according to Eqs. IjlGI) and 
()19|l - i.e., it grows during the expansion - which is (al- 
most) certainly not a realistic feature of the real Planck 
length. In order to circumvent this obstacle, one could 
introduce another cut-off scale which is sensitive to fre- 
quencies instead of wavenumbers. For example, via ex- 
posing the condensate to monochromatic radiation with 
a detuned radio-frequency lu = loq — Aco with luq being 
the frequency of a suitable hyperfinc transition and Auj 
the detuning, one may couple frequency-sclcctivcly cou- 
ple phonons with the frequency Alo to that transition. 
In that case, Auj would be a reasonable analogue of the 
Planck scale, by fixing a frequency scale at which free 
phonon propagation is cut off. 



D. Outlook 

After succeeding to measure the main effect - the fact 
that quantum fluctuations of the inflaton fleld in de Sit- 
ter space-time get frozen at a specific value - we can 
start to manipulate the behavior at the effective Planck 
scale and study its impact on the frozen fluctuations. 
The exciting prospect of undertaking such manipulations 
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of Planck scale physics is that one were able to experi- 
mentally investigate analogue signatures of new (trans- 
Planckian) physics in the anisotropies of the cosmic mi- 
crowave background; such possible signatures of trans- 
Planckian physics are discussed, e.g., in 31j. 

We, finally, mention that a further interesting aspect 
of the scale separation of independently propagating spin 
and density modes is the possible implementation of a 
variant of the supersonical "warp-drives" proposed in 
[33 • The reference frame, against which the "superlumi- 
nal" motion of the 0_ mode is measured, is not provided 
by the laboratory frame, like in |32|. but by the (/>_)_ mode, 
which has effectively instantaneous signal transfer on the 
time scales of the </)_ mode. In addition, the "hard" 0+ 
mode is essentially flat on the curvature radius scale of 
the (/)_ mode, i.e., appears to be essentially Minkowski 
space-time as seen from the 0„ mode. The Minkowski 
background of the 0+ mode represents a necessary pre- 
requisite for "superluminal" travel in warped space-times 



to be operational and definable, because it enables the 
comparison of two metrics, one of them flat and the 
other curved, on the same manifold [23 • In contrast 
to the single Euler- fluid case of '33| , in the presently dis- 
cussed two-component Bose-Einstein-condensed gas both 
of these metrics obey Lorentz invariance. 
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